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ABSTRACT 

Direct observations of the supposedly universal primordial deuterium 
abundance imply a relatively large baryon density Qb = (0.019 — 0.030)/i“^ 
(95% C.L.). On the other hand, concordance between the previously accepted 
and ^Li abundances and standard Big Bang Nucleosynthesis requires the 
thrice smaller value = (0.005 — 0.010)/i“^ (95% C.L.). For each Qb, we use 
X-ray and Sunyaev-Zeldovich observations of the baryon fraction Jb in rich 
clusters of galaxies, in order to obtain limits on the total mass density Qci in 
clusters of galaxies. The higher-De values are consistent with clusters being 
a fair sample of the Universe, and then imply = (0.3 — 0.9), a medium or 
critical density Universe. Said otherwise, the observed limits /b > 0.1, Qm > 0.3 
imply Qb > 0.03. If the newer^Re abundance observations are accepted, this is 
consistent with standard BBN. 
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1. DIFFERENT BARYON/PHOTON RATIOS OBTAINED FROM BBN 


1.1. The Former Crisis in Standard Big Bang Nucleosynthesis 


Deuterium has been traditionally used as a baryometer. Until recently, observations 
of the chemically evolved nearby interstellar medium and solar system determined only a 
lower bound D/H > (1.6 ± 0.2) x 10“^. In the standard three low-mass neutrino scenario 
for Big Bang Nucleosynthesis (SBBN), this determines an upper bound tothe primordial 
mass ratio, t]io = < (8.2 ± 0.5) (|Hata et a/.(1997)|, [Bludman (1996)|) . Recent 

observations ( P?ytler et a/. (199^ [Buries and Tytler (1996)|) of deuterium absorption lines 
in two nearly primordial Lyman limit clouds illuminated by distant quasars (QSA) show 
D/H = (2.4 ± 0.3 ± 0.3) x lO"® or r^io = 6.4tE];?. 

Both these high values for r]iQ are inconsistent with the low value r^io = 1.8 ± 0.3 once 
reported ([Carswell et al. (1994)| , [Songaila et a/.(1994)[ , [Rugers and Hogan (1996)1) in light 
from quasar Q0014-I-813, but since discredited ([Tytler et a/.(1997][) . This low rjio had been 
in excellent concordance with the abundance. 


Yp = 0.232 ± 0.003(stot) ± 0m5{syst) 


( 1 ) 


( [Geiss (1993)[, Plive and Steigman (199^) , formerly inferred from HII regions and the 
Standard = 3 model) and in good agreement with the Spite plateau in halo 
stars (^pite et a/. (19841) . The higher rjio value, on the other hand, was inconsistent with 
the above ^He abundance and demanded depletion in stars and surprisingly little 
Galactic chemical evolution of ^He. This discordance between the higher rjiQ = 6.41°; 
demanded by low D/H and the lower r/io = 1.8 ± 0.3 demanded by the then-accepted 
Helium abundance (1) was the crisis for BBN ([Hata et al. (199^). 


After this paper was hrst submitted for publication, improved measurements of D/H 
and Yp have ruled out the lower t]iq value and ended the BBN crisis: (1) Improvements 
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in the two QSA determinations (|Tytler et a/.(1997]|) now give D/H = (3.2 ± 0.4) x 10“^, 
implying = 0.020 ± 0.002, r^io = 5.5 ± 0.5, Yp = 0.247 ± 0.0025; (2) This is consistent 
with Yp = 0.243 ± 0.003 measured ( [Izotov and Thuan (19971 ) in an important Wolf-Rayet 
galaxy I ZW 18 and in eight other low-metallicity HII clouds. These increases m D/H and 
in Yp suggest concordance at values for t^iq, and flci only 15% lower than the 

values we adopted in the last column of Table 1, sustaining our conclusions in Section 4. 


1.2. Implications for the Baryon and Total Mass Densities in the Universe 


Because the photon number density in the Universe is well determined by the cosmic 
background temperature, the two incompatible values t]io = or (1.8 ±0.3), determine 

substantially different baryon mass densities Qb = 0.0036677ioh“^, where h ~ 0.70 ± 0.15 is 
the present Hubble constant in units of 100 km/sec/Mpc ([Freedman et al. (1994)| , [Riess[ 
et a/.(1995i [Kennicutt et al. (1995)[) . For the allowed range Hq = (70 ± 15) km/s/Mpc, 


either Qb = (0.005 — 0.010) or (0.019 — 0.030) (95% C.L.) ([Hata et a/.(1997)|) , which we 
hereafter we refer to as lower- and higher-Hs (Table 1). Unless h is large, the lower-H^ 
may already be excluded by present baryon inventories ( [Fukugita et al. (1996)[ , [Bahcall 


[( 1997)[) showing > 0.03. In the next section, we discuss two different measures of the 
baryon mass fraction /s = Mb/Mci in rich clusters of galaxies and use them to calculate 
the total mass density in rich clusters, Qci = ^b/Jb- In Section III, we calculate the limits 
on Qb/Jb for the lower- and higher-Hs separately. Because Jb in rich clusters might be 
enhanced over the cosmic value Qb/^tti by a “baryon enhancement factor” T ([White ej 
a/.(1993)[, [Steigman and Felten (1995[), the cosmic total mass density flrn = ^^b/Ib might 


differ from Qci defined above. 


This paper extends Hata et al. (1997) by including new data ([Evrard (1997)[ , [Herbig elj 
al. (1995[, [Myers et al. (1997)[) on the baryon fraction in rich clusters and new dynamical 
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limits on et aL (1997)|) , and concludes that clusters are indeed a fair sample of 

the Universe (T ~ 1), that the baryon density is relatively high {Qb = (0.025 — 0.099) (95% 
C.L.), that h ^ 0.55 allows a critical density Universe, but that larger h implies that we live 
in a medium density Universe. 


2. BARYON MASS FRACTION IN RICH CLUSTERS 


The luminous matter density Vtium = 0.004+ 0.007h ([Fukugita et al. (1996)| , [Bahcall 
EMIE ersic and Salucci (1992)|) . The observed mean Lya flux decrement shows baryon 
density VLb > 0.0125h“^ ([Weinberg et al. (1997)| . This shows that considerable baryonic 
matter is dark and that SBBN predicts D/H < 6 x 10“®. This already argues against the 
lower rjiQ or choice 


In this section, we will summarize dynamical observations showing that, the baryon 
fraction in rich clusters, fs < 0.18, so that rich clusters are dominated by non-baryonic 
matter. These rich clusters are the largest virialized structures and, although of intermediate 
size ((1 — 10)h“^ Mpc) and constituting only a small fraction of the total mass in the 
Universe, are thought to be fair samples of the entire Universe. 


2.1. Baryon Fraction in Cluster Hot Gas 


The baryonic and total masses, M^, Md in the hot gas in clusters have each been 
measured by two different methods: originally ([White and Frenk (1991)| , [White and Fabicin 
(1995 )[, Evrard (1997)[) from the X-ray bremmstrahlung off hot cluster gas, and, more 


recently ([Myers et al. (1997)[), from the Sunyaev-Zeldovich (SZ) inverse Compton spectral 


distortion of cosmic background radiation. Both methods depend on modeling the density 
and thermal structure of clusters to determine the baryon mass in gas Mb, the total mass 
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Mci, and, hence, the clnster baryon mass fraction fs = Mb/M^i = Qb/^cI- 


From detailed analysis of a large sample of rich clnsters, and comparison with CDM 
clnster models, Evrard ([Fvrard (1997)) obtains for the baryon fraction in hot gas 


fB{Xray, Evrard) = (0.060 ± 0.003) h 


( 2 ) 


This is consistent with 


fB{Xray, select) = (0.054 ± 0.013) h 


(3) 


observed in three large, well-stndied clnsters A2142 ([White and Fabian (1995)1) , A2256 
( [Henry et al. (1993)|) and Coma ( [White et a/. (199^) , that are probably ronnder and 
smoother than the other clusters in Evrard’s larger, less selective catalogue. (Both these 
X-ray measurements are consistent with earlier, less sensitive measurements ([White and 
Frenk (1991)[) 

/B(Xra|/) = 0.049l°;°28h-3/2, (4) 


which we do not show or use.) 


From the Sunyaev-Zeldovich effect in these three rich clusters, Myers et al ([Myers ei 
al. (19971 ) obtain 

Jb^SZ) = (0.061 ± 0.011) h-\ (5) 

Both values (2) and (4) are unweighted means, with the errors combined in quadrature. 
(Because a fourth cluster, A478, has a SZ mass three times higher than these three and 
differs signihcantly in other respects, we and Myers omit it from our determination of 
Jb- Had Myers included A478, he would have obtained fB{SZ) = (0.087 ± 0.030) h~^, 
insignihcantly higher than the baryon fractions (2) and (4). Besides A478, we have also 
omitted from the SZ analysis, clusters A665 ([Birkenshaw et al. (1991)[) and CL0016-I-16, 
A773 ([Carlstrom et al. (1996)[) for which X-ray observations are apparently not available.) 
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2.2. Consistency of X-ray and SZ Determinations of Baryon Fraction 

For these three selected clusters, Myers et al obtain SZ hot gas masses that are on 
average (1.27 ± 0.13)h^/^ larger than the X-ray hot gas masses measured by White and 
Fabian ([White and Fabian (1995)|) , by Henry et al ([Henry et al. (199^ , and by White 
et al ([White et a/. (199^ ), i.e. the X-ray and SZ hot gas masses from these three nearby 
clusters are consistent at 

fsiHG, select) = ( 0.12 ± 0 . 02 ) ( 6 ) 

over the broad range h = 0 . 62 lg;]^ 2 - smaller h, this selected cluster baryon fraction 
is some 11% less than that observed in Evrard’s larger X-ray sample. Although Myers 
excludes A478 from his determination of h, when he includes it, he obtains h = 0.54]'^q;J^, 
lower than but consistent with the value we obtained from the three selected clusters. 


2.3. Baryon Fraction in Galaxies 


The baryons in hot gas must be augmented by the luminous mass flium ~ 0.009 
( [Lynds and Petrosian (1986)[ , ^oucail et a/.(198'7][ , [Tyson et a/.(1990][, [Fukugita et al. 


[(1996)[, [Bahcall (1997)[) and by dark baryons residing in galaxies. Judging by the halo 
mass in our own Galaxy, the dark mass fraction could be as large as this luminous mass 
fraction, making Mgai/Md = 0.0135 ± 0.0045. (Myers adopts the slightly smaller value 
Mgai/Mci = 0.009 ± 0.003, observed in Coma ( [White et a/.(1993')[) ). We therefore hnallly 
adopt (Fig. 1) for the total baryon mass fraction, as measured in X-rays, 


feiEvrard) = (0.060 ± 0.003) + 0.0135 ± 0.0045, 


(7) 


or 


fBiselect) = (0.054 ± 0.013) ^ o.0135 ± 0.0045 


( 8 ) 





























and, as measured by the Sunyaev-Zeldovich effect, 


fsiSZ) = (0.061 ± 0.011) h-^ + 0.0135 ± 0.0045. (9) 

The Icr bounds on these three cluster total baryon mass fractions are shown in Figure 1 by 
solid and by dashed lines for Evrard’s and for the three selected X-ray clusters respectively, 
and by the shaded area for the SZ observations. For our preferred range h = 0.621qj2, 

Jb ~ 0.13 ± 0.03. For larger h, the three observations agree on fb ~ 0.09 ± 0.02. For 
smaller h, the larger values fb ~ 0.12 ± 0.02 from SZ and ft, ~ 0.15 ± 0.03 from X-rays 
are suggested. Note that including the X-ray measurements and the baryons in galaxies 
increases the baryon fraction over that in Eqn. (4). In any case, the conservative bounds 
are fs ~ 0.11 ± 0.03 for larger h and fs ~ 0.14 ± 0.04 for smaller h. 

These upper and lower bounds on the baryon fraction provide lower and upper 
limits on fls/fB = ^ci, shown in Figure 2 for lower-(solid curves) and for higher- 
(dashed curves). The lower = (0.007 ± 0.001), implies Qci < 0.27. The higher 
= (0.023 ± 0.003) implies = (0.3 — 0.9). This estimate is a little smaller than 
the value obtained ([Myers et al. (1997]|) using SZ data alone. Nevertheless, is still 

possible, for h < 0.6. 

We recall that baryon inventories (jFukugita et al. (1996)| , Pahcall (1997)|) already 
show > 0.03 and that the cosmic virial theorem already implies Qm ^ 0.2. If we were 
ready to accept rich clusters as a fair sample of the Universe, this would already practically 
exclude the lower value. Nevertheless, in the next section, we will test the possibility 

that the cluster and cosmic total mass densities differ by a baryon enhancement factor 
T = 
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3. DYNAMICAL MEASURES OF THE COSMIC MASS DENSITY 

If T ~ 1, the lower and higher VLb choices discussed in Section II would then already 
imply respectively a low and an intermediate or critical density Universe. Although 
numerical simulations and other theoretical arguments ([White and Frenk (199lj| [White 
[and Fabian (199^ strongly suggest that rich clusters do not appreciably concentrate 
baryons, we will now test this assumption by comparing VLd from clusters with large-scale 
determinations of Most of these large-scale determinations depend on models for the 
evolution of large scale structure from assumed initial fluctuations and dark matter content 


The least model-dependent global bounds on derive from (1) diverging flows in voids, 
(2) from distant {z ~ 0.4) supernovae la distance indicators, (3) from weak gravitational 
lensing of quasars by intervening galaxies and rich clusters, and (4) from the expansion age 
of the Universe, to = Ho^ From diverging flows, ilm > 0.3 (2.4ct) because “voids 

cannot be more empty than empty” (mass densities cannot be negative) (Pekel and Rees[ 
(1993)[) . From type la supernovae (Perlmutter et q/.(1997)[) , —0.3 < < 2.5 (90% 


C.L.) or Qrn > 0.49 (95% C.L.) for a flat Universe. From the mass-to-light ratios in weakly 
lensed rich clusters ([Kaiser (1995)[) , Qm — (0.3 — 1) and from the statistics of gravitational 
lens counts ( Kochanek (19961) , > 0.34 in a flat universe From the most likely age of the 

oldest stars to > 12 Gyr ( [Ghaboyer et al. (1997)| , [Keid (1997)| ), implies the upper bounds 
on Gm shown by the heavy dashed curves in Figure 2, for an open cosmology with Ga = 0 
and for a flat cosmology with Ga = 1 ~ If T = 1 and to, h, fs are all at their lower 
limits, a matter-closed universe Gm = 1 is possible. 


On large comoving scales (10 — 100)h“^Mpc, the total mass density derives from 
galaxy redshift surveys (subject to optical biasing), from dynamical studies of cosmic flows, 
and from GBR growth of fluctuations (free of optical biasing) . Omitting the observations 
which depend on optical biasing, Dekel et al (Pekel et al. (1997)[) summarize: (1) From 
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preliminary CAT and Saskatoon observations ([Hancock et at. (1996)|, [Nettertield et at. 


|( 1997)1) of the first CBR acoustic peak, + Ha > 0.3 (95% C.L.); (2) Cosmic Flows 
( ^aroubi et a/.(1997)|, [Kolatt et al (1997)|) give > 0.3 ; (3) Growth of Fluctuations 
together with cluster morphology give Qm > 0.2; (4 ) Cobe Power Spectrum + Mark III 
Velocities gives model-dependent results: (a) ~ 0.2h~^ from untilted CDM models, 

assuming no optical biasing; (b) ~ (0.45 ± 0.07) from spatially flat CDM with spectrum 

tilt n; (c) the best CDM £t requires a small tilt and either ~ 0.7 or ~ 0.2; (d) the 
hrst acoustic peak requires a small tilt and a high baryon content, Qb ~ 0.1 ( Zaroubi~ 
a/.( 1997)1) . 


All these observations of large-scale structure require > 0.3, so that since 

Jb = 0.14 ± 0.04 (0.11 ± 0.03) for small (large) h, we have > 0.03 (0.024). Unless rich 

clusters are baryon-enhanced by the unreasonable factor T ~ 3, this rules out the lower-Hg 
solution, Hs < 0.01 (95% C.L.). The large-scale structure observations are consistent with 
the higher-Hs observations, = (0.025 — 0.099) (95% C.L.) and with rich clusters being a 
fair sample of the Universe. If h is near its lower bound, and VLb near its upper bound, a 
critical density universe = 1 is just allowed. 


4. CONCLUSIONS 

We now summarize (Table 1) our cosmological conclusions, distinguishing the 
implications of the baryon fraction observed in X-ray and SZ studies of rich clusters, from 
the implications of large-scale structure. In all cases, we assume Hq = (55 — 85) km/s/Mpc. 

The lower value f^B = (0.008 — 0.033) (95% C.L.) together with the baryon fraction in 
rich clusters would have implied a low cluster mass density Qci < 0.27, inconsistent with 
baryon inventories and with rich clusters being a fair sample of the Universe. 
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The larger value = (0.025 — 0.099), together with the baryon fraction observed in 
both X-ray and SZ clusters, allows a higher matter density Qci = (0.3 — 0.9), consistent with 
rich clusters being a fair sample of the Universe and with data from large-scale structure. 
If Us is near its upper bound and h near its lower bonnd, a critical density is possible. 

This work is snpported by Department of Energy Contract No. DE-AC02-76-ERO-3071 
and has benehted by spirited discussions with G. Steigman. 
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Table 1: Cosmological implications of lower-(consistent with previonsly observed 

'^Li abnndances) and of higher-fl^ (consistent with the denterinm abnndances observed 
in the nearby Galaxy and in two qnasar aborption systems ([Tytler et a/. (1996)], Pnrles and 
[lytler (1996)1 , |Geiss (1993)| , |Hata, bcherrer et al. (199F)|) . The errors are for 68% C.L., 
while the ranges in the parentheses are for 95% C.L. We take Hq = (70 ± 15) km/s/Mpc, 
so that = (0.49 ± 0.21). In each case, the mass density in clnsters, kid = kls/fs, is 
obtained from the range in fs in Fignre 1. (If the latest determinations ([lytler et al. [1997) ) 
D/H = (3.2 ± 0.4)10“^, r]io = 5.5 ±0.5 were now nsed, kid in the last colnmn and 

in Fig. 2 wonld now be rednced by 15%.) 



lower-D b 

higher-D b 

D/H (lO-*^) 

19 ±4 

2.4 ±0.3 ±0.3, > 1.6 ±0.2 

^10 

1.8 ±0.3 (1.7- 2.7) 

6.4lg;? (5.1 - 8.2) 

kl^h^ 

0.007 ±0.001 

0.023 ±0.003 


(0.005 - 0.010) 

(0.019 - 0.030) 

kls 

(0.008 -0.033) 

(0.025 -0.099) 


klcl 


0.09 -0.26 


0.3 - 0.9 























17 


Fig. 1.— The allowed la range of baryon fraction fs in rich clnsters of galaxies as fnnction 
of the Hnbble constant Hq = lOOh km/s/Mpc, as measnred by thermal bremmstrahlnng 
X-rays (in Evrard’s catalogne( |Evrard (1997)|) (thin solid cnrve) and in three selected rich 
clnsters ( [White and Fabian (1995)| , [Henry et at. (199¥)| , [White et a/.( 1993]]) (dashed cnrve)) 
and by the Snnyaev-Zeldovich npscattering of cosmic backgronnd radiation (shaded region 
SZ) ([Myers et al. (1997)[) . In all three cases, the baryon fraction measnred in hot gas has 
been angmented by the Inminons and dark baryons in galaxies and stars. 


Fig. 2.— BBN and clnster baryon fraction constraints on the clustered matter density 
VLci = I^b/Jb as function of Hubble constant h for = 0.007 ± 0.001 (solid curves) and 
for = 0.023 ± 0.003 (dashed curves). The regions between each pair of curves are 

allowed at 68%(95%) C.L. The heavy dashed curves are upper bounds on the cosmic baryon 
density in open and flat Universes of age greater than 12 Gyr. In principle, niight 
differ from the clustered baryon density Qci by a factor T = Qm/^ci- In fact, observation of 
large-scale structure show Qm ~ ^ci, so that T ~ 1: rich clusters are a fair sample of the 
Universe. 
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